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ABSTRACT 

TEMPERATURE  AND  PRESSURE  COMPONENTS 
OF  SHOCK-INDUCED  CHANGES  IN  THE 
ELECTRONIC  SPECTRUM  OF  CARBON  DISULFIDE 

by  Richard  H.  Granholm,  M.S. 

Washington  State  University,  August  1985 

Chairman:  G.  E.  Duvall 

Liquid  carbon  disulfide  has  its  first  electronic  transition  centered  at  3200A  . 
Shock-compression  of  CS2  up  to  about  80  kilobars  has  been  shown  to  shift  the  red  edge 
of  this  absorption  band  more  than  1600 A  towards  longer  wavelengths.  In  this  work  the 
conditions  of  temperature  and  pressure  behind  the  shock  wave  are  controlled  to  allow 
their  separate  effects  on  the  CS2  spectrum  to  be  determined.  The  shift  of  the  band  edge 
appears  to  be  a  function  only  of  temperature  and  pressure. 
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SECTION  1 

INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 

It  has  been  known  since  1981  that  liquid  CS2  exhibits  a  dramatic  shift  in  the  red 

O  1 

edge  of  its  3200A  absorption  band  when  subjected  to  shock-wave  compression. 

Perhaps  the  earliest  chemical  reaction  initiated  by  shock  waves  was  the  formation  of 
organic  molecules  in  the  origin  of  life  on  Earth.  More  recently,  interest  in  the  chemical 
effects  of  shock  waves  has  been  fueled  by  reports  from  Russian  scientists  on  unusual  pro¬ 
ducts  recovered  from  shocked  materials,^’^  and  by  technical  advances  in  this  country 
on  dynamic  measurements  in  the  time  scale  of  shock  experiments.  These  new  techniques 
make  it  possible  to  study  the  early  steps  in  shock-induced  reactions. 

Bridgman  first  reported  a  reaction  in  liquid  CS2  when  it  is  compressed  statically  to 

>45  kbar  (1  kbar  *  987  atm.)  and  175°C.^  Later  work  suggested  the  formation  of  a 

7  8 

carbon-sulfur  polymer  from  CS2  under  similar  conditions.  ’  With  velocity  measurements, 
Dick  found  a  transition  to  a  higher  density  state  in  CS2  occurring  in  shock  experiments  at 
about  60  kbar.^’1®  Sheffield11  studied  the  kinetics  of  this  shock-induced  transition,  and 
found  an  induction  period  preceding  relaxation  to  the  new  state,  which  suggested  a  mul¬ 
tistep  mechanism. 

The  first  spectral  measurements  on  CS2  at  high  pressures  were  made  by  Ogilvie  and 
Duvall.1  Their  first  experiments  showed  total  extinction  of  light  transmitted  through  a 

O 

bandpass  from  2670  to  2840A  between  two  absorption  regions  in  the  CS2  spectrum  (see 
figure  1).  The  shock  pressures  ranged  from  83  to  127  kbar.  They  later  found  that  a  shock 
of  about  5-10  kbar  was  sufficient  for  extinction  of  this  bandpass;  but  more  significantly, 
the  extinction  occurred  by  apparent  shifting  or  broadening,  with  pressure,  of  the  lower 
wavelength  absorption  band  towards  the  red.  The  320oA  absorption  band  also  expanded 
towards  the  red  with  increasing  pressure. 

The  goals  of  the  present  work  are  to  measure  this  shock-induced  shift  in  the  CS2 
absorption  spectrum,  and  to  sort  out  the  effect  of  temperature  from  that  of  pressure. 
These  goals  were  met  by  devising  experiments  in  which  the  sample  temperature  behind 


the  shock  wave  can  be  modified,  thus  allowing  the  wavelength-P-T  surface  for  CS2  to  be 
mapped  out  This  mapping,  i.e.,  the  results  of  the  experiments  in  this  report  is  shown  in 
figure  2. 

The  present  work  deals  only  with  "unreacted"  CS2,  that  is,  with  changes  that  occur 
prior  to  reaching  the  temperature  and  pressure  threshold  of  the  transition  observed  by 
Dick  and  by  Sheffield. 

1.2  SUMMARY 

Liquid  Carbon  Disulfide  exhibits  a  dramatic  shift  in  the  red  edge  of  its  3200A 
absorption  band  when  compressed  by  shock  waves.  In  the  present  work  this  shift  was 
studied  at  various  temperatures  and  pressures  from  201°AT  to  ~750°K  and  from  1  atmo¬ 
sphere  to  ~90  kilobars.  These  conditions  were  reached  in  shock  experiments  with  heated 
and  cooled  spectroscopy  cells,  so  that  the  CS2  spectrum  could  be  observed  at  a  given 
pressure  with  different  temperatures,  or  at  a  given  temperature  with  different  pressures. 
From  analysis  of  these  experiments,  the  spectral  shift  appears  to  be  a  function  only  of 
temperature  and  pressure.  These  two  factors  appear  to  be  independent  dX/dT  -vs-  T  is 
approximately  linear,  but  dX/dP  shows  no  clear  dependence  on  either  temperature  or 
pressure,  but  appears  to  remain  constant  within  ~±  20%. 


Figure  2:  P-T-A  Surface  For  Liquid  CSj 

CS  IhiaS  SJ?wsJ^e  positjon  of  the  red  ed8«  of  3200A  absorption  band  in 

i  pTifSure  and  temperature.  The  three  curves  depicting  the  surface 

201o?  293°J?  and  lolS  00  CS1  at  three  djfferent  initial  temperatures: 

ui  a.,  2V3  K,  and  398  K.  The  dotted  lines  connect  points  of  constant  pressure. 
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2.1  MAIN  FACILITY 

The  experiments  were  conducted  with  the  4  inch  gas  gun  at  the  Shock  Dynamics 
Laboratory,  Washington  State  University.  The  gun  is  46  feet  long  and  has  a  compressed- 
gas  breech  useable  to  6000  psi.  Velocities  up  to  about  1.2  km/sec  can  be  attained  with  a  1 
kilogram  projectile.  The  target  chamber  contains  optical  and  electrical  ports  for  recording 
the  experiments,  and  a  rag-filled  catcher  tank,  contained  within  a  concrete  bunker.  A 
detailed  description  of  this  facility  can  be  found  in  reference  14. 

A  gas  gun  is  well  suited  for  high-pressure  research  since  it  can  generate  a  planar 
shock  wave  that  is  well-characterized,  allowing  the  pressure  history  of  a  sample  to  be 
accurately  measured  or  calculated. 

2.2  SPECTROSCOPIC  APPARATUS 

A  block  diagram  of  a  spectroscopy  experiment  is  shown  in  figure  3.  The 
spectrograph-streak  camera  consists  of  a  Beckman  &  Whitley  model  1500  grating  disper¬ 
sion  unit  and  model  339B  continuous  writing  streak  camera  using  35mm  film.  A  diagram 
of  the  recording  system  is  shown  in  figure  4.  A  detailed  description  can  be  found  in 
reference  15. 

Figure  5  shows  a  sketch  of  the  spectroscopy  cell  target  assembly  and  a  projectile. 
Sapphire  is  used  as  the  window  material  because  it  remains  transparent  at  the  pressures 
reached  in  these  experiments.  The  impactor  and  back  window  are  one  inch  in  diameter 
and  ~l/2  inch  thick.  The  front  window  is  1  and  1/4  inches  in  diameter  and  typically  1/8 
inch  thick.  The  sample  layer  is  typically  150  microns  thick.  The  cell  is  filled  through 
1/16  inch  copper  tubes.  The  light  source  is  a  xenon  fiashlamp  which  is  mounted  on  the 
face  of  the  target  and  is  flashed  when  a  trigger  pin  on  the  target  is  shorted  by  the  projec¬ 
tile.  The  entire  assembly  at  the  moment  of  impact  is  shown  in  cross  section  in  figure  6. 
Projectile  and  target  design  are  by  Paul  Bellamy^. 


SPECTROGRAPH 


Figure  3:  Outline  Of  A  Spectroscopy  Experiment 

The  xenon  flashlamp  is  triggered  just  prior  to  impact  The  light  passes  through  the  sap¬ 
phire  impactor,  A,  then  through  the  sample  cell  mounted  in  the  target  B  and  C  are  the 
front  and  back  sapphire  windows  of  the  sample  ceiL  The  light  is  then  directed  to  a  spec¬ 
trograph  where  a  time-resolved  spectrum  is  recorded  as  the  sample  is  compressed  by  the 
shock  wave. 


sapphire  front  window;  B)  sapphire  back  window.  (Drawing  by  S.  Ar 


Figure  6:  Projectile  And  Target  Assembly  At  Impact 
(Design  and  drawing  by  P.  Bellamy.) 


SECTION  3 

DESCRIPTION  OF  EXPERIMENTS 
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3.1  SHOCK-COMPRESSION  OF  CS2  :  STEPWISE  LOADING  (SWL) 

When  the  impactor  strikes  the  front  window  of  the  spectroscopy  cell,  a  shock  wave 
is  propagated  forward  into  the  window  (and  a  backward-facing  wave  propagates  into  the 
impactor),  causing  a  sudden  jump  in  pressure  and  temperature.  When  this  wave  enters  the 
liquid,  its  amplitude  is  reduced  because  of  the  lower  shock  impedance  of  the  liquid  com¬ 
pared  to  sapphire,  but  when  the  wave  reaches  the  rear  window  it  is  reflected,  causing  a 
second  jump  in  pressure  and  temperature.  In  this  way,  the  pressure  in  the  liquid  "rings 
up"  through  a  series  of  reflections  until  it  reaches  the  initial  sapphire-sapphire  impact 
pressure.  This  process,  which  we  shall  call  stepwise  loading  (SWL),  is  illustrated  in  pres¬ 
sure  -  particle-velocity  space  in  figure  7. 

The  normal  spectrum  of  liquid  CS2  is  shown  in  figure  1,  and  a  time-resolved  spec¬ 
trum  of  CS2  in  a  shock  experiment  is  shown  in  plate  3.  Xenon  lines  are  visible  on  the  left 
in  plate  3,  until  the  gas  in  the  flashlamp  becomes  hot  enough  for  broad  band  emission. 
The  3200A  absorption  band  appears  at  the  bottom  of  the  record,  with  its  upper  edge  at 
3550A .  At  impact  there  is  a  slight  increase  in  total  transmitted  light  when  the  reflective 
surfaces  of  impactor  and  cell  window  are  merged.  280  nanoseconds  later  the  shock  enters 
the  CS2  causing  a  sudden  jump  to  8.5  kbar  and  428° AT.  Simultaneously  the  edge  of  the 
absorption  band  shifts  about  200 A  to  the  red.  The  shock  wave  traverses  the  layer  of  CS2 
(typically  150  microns)  in  ~85ns  at  which  point  it  is  reflected  back  into  the  CS2  causing 
another  jump  in  temperature  and  pressure  and  another  shift  of  the  band  edge.  The  first 
three  or  four  such  steps  are  visible  in  this  record.  When  the  final  values  are  reached,  81 
kbar  and  700°X  ,  the  band  edge  levels  off  at  ~5300A . 

Pressure  and  temperature  are  difficult  to  measure  directly  in  shock  experiments  so 
these  values  are  calculated  for  each  step  from  an  equation  of  state  (see  section  4).  Even 
after  the  steps  become  too  small  to  be  seen  on  the  record,  their  location  in  time  is  known 
from  the  known  values  of  shock  velocity  in  CS2-  Thus  one  experiment  can  give  10  or  12 
points  on  a  curve  of  wavelength  -vs-  pressure.  The  Pressure  history  of  the  sample  in  a 
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PRESSURE 
k  bars 


mm/^s 


Figure  7:  Stepwise  Loading 

The  shocking  up  process  in  a  CS2  layer  between  sapphire  plates.  The  slope  of  each  curve 
is  a  characteristic  of  the  materiar(sapphire  or  CS2).  Point  A  is  the  initial  projectile  velo¬ 
city,  and  0  is  the  initial  target  condition.  Point  B^represents  the  conditions  of  the  initial 
sapphire  on  sapphire  impact.  The  reverberating  shocks  take  the  CS2  to  states  1,  2,  3,  etc. 


typical  experiment  is  graphed  in  figure  8.  Alignment  of  target  and  projectile  is  critical  in 
these  experiments  since  tilt  would  introduce  error  in  the  calculated  position  in  time  of  the 
shock  wave,  by  "smearing  out"  the  stepwise  jumps  in  pressure,  temperature,  and 
wavelength. 

3.1.1  ROOM  TEMPERATURE  SWL  EXPERIMENTS 

In  all  the  shock  experiments  described  in  this  report,  the  original  experimental 

record  is  a  strip  of  35mm  film  with  the  CS2  spectrum  dispersed  in  wavelength  (A)  across 

the  film  and  in  time  along  the  film,  as  seen  in  plates  1  through  7.  The  first  step  in  analysis 

0 

is  digitizing  the  position  of  the  3200  A  band  edge.  This  is  accomplished  on  a  Gaertner 
measuring  microscope  interfaced  to  an  Apple  computer  for  recording  the  data.  A  fortran 
program  is  run  separately  to  simulate  the  experiment  This  program  is  called  SHOCKUP 
and  was  originally  written  by  Kendal  Ogilvie13  to  calculate  pressure,  temperature,  and 
time,  using  the  ShefTield-Duvall  equation  of  state  (see  section  4)  and  the  shock  jump  con¬ 
ditions,  as  the  CS2  "shocks  up"  in  a  SWL  experiment  The  program  currently  in  use  was 

1  7 

written  by  Robin  Collins. 

The  original  data  give  A  -vs-  time  and  SHOCKUP  gives  P  -vs-  time.  By  matching 
the  time  in  the  experiment  to  the  time  calculated  with  the  program,  we  obtain  AA  -vs-  P, 
where  AA  is  the  shift  of  the  band  edge  relative  to  its  original  position.  The  original  shot 
record  of  plate  3  is  shown  in  digitized  form  in  figure  9.  The  times  of  each  stepwise  jump 
in  A  are  readily  matched  with  SHOCKUP  times,  which  are  plotted  on  the  same  graph. 

Five  room-temperature  experiments  were  performed,  i.e.,  the  CS2  samples  were  ini¬ 
tially  at  ~20°C.  The  impact  pressures  ranged  from  54  to  130  kbars.  (One  of  these  experi¬ 
ments,  shot  82-015,  was  performed  by  Paul  Bellamy;  the  rest  were  by  this  author.)  Copies 
of  the  original  shot  records  are  shown  in  plates  1  through  5.  The  AA  -vs-  P  data  from  all 
five  experiments  are  plotted  in  figure  10,  and  are  fit  quite  well  by  the  cubic  equation: 

AA  -  26.75P  -  0.2554P2  +  0.002252P3  (AA  in  A;  P  in  kbar) 


as  shown  by  the  solid  line  in  the  figure.  The  T-P  profile,  calculated  with  SHOCKUP,  for 


PRESSURE  (kbar) 


TIME  (nanoseconds) 


Figure  8:  Pressure  History  Of  CS2  In  A  Typical  Experiment 

This  shows  the  pressure,  just  behind  the  reverberating  shock  wave,  of  a  160  micron  thick 
sample  of  CS*  in  a  stepwise  loading  experiment.  The  shock  just  enters  the  liquid  at  time 
*  0.  Impact  pressure  is  8 1  kilobars. 


DELTA  LAMBDA  (A) 
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TIME  (nanoseconds) 


Figure  9:  Shot  82-037;  TQ  *  20°C  (293° K ) 

O  -  AA  -vs-  time;  digitized  band  edge  from  original  shot  record. 

X  -  Pressure  -vs-  time;  calculated  by  SHOCKUP. 

"  0?  the  shock  just  enters  the  CS*.  At  the  first  X,  the  shock  completes  one  transit 
ot  the  CS-j  and  is  just  about  to  reflect  frdm  the  back  sapphire  window.  At  this  instant  the 
”2  Prcs.  w  homogeneous  throughout  the  cell.  The  dotted  line  shows  how  pressure, 
calculated  by  SHOCKUP,  is  correlated  with  the  wavelength  shift  of  the  band  edge. 
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Figure  10:  AA  -vs-  P  For  Room  Temperature  SWL  Experiments  (Tq  *  20°C) 

The  circles  show  the  position  of  the  3200A  CS^  band  edge  for  each  pressure  step  in  the 
five  shock  experiments  recorded  in  plates  1-5.  The  solid  line  is  a  least  squares  fit  of  the 
data.  The  equation  of  the  line  is: 

AA  -  26.75P  -  0.2554P3  +  0.002252 P3 

AA  is  in  angstroms  and  P  is  in  kilobars. 
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Figure  11:  T-P  Profiles  Calculated  By  SHOCKUP 

Curve  A)  cooled  cell  experiment,  7*o  ■  201°K;  curve  B)  room  temperature  experi¬ 
ments,  T0  *  293 °K\  curve  C)  heated  cell  experiment,  T0  *  398 °K.  The  circles  are  the 
actual  SHOCKUP  data;  the  solid  lines  are  cubic  fits  to  the  data.  The  equations  for  the 
(inverse)  lines  are  as  follows: 

A)  P  -  -8.075  +  5.358x  10~3r  -  2.429x  lO^T3  +  8.765x  10~77^ 

B)  P  -  -32.94  +  2.544x  lO^T  -  7.378x  [Q^T2  +  8.640x  10"77^ 

C)  P  - -185.6  +  1.2007- 2.61  lxlO^T2  +  1.930x10^ 

where  P  is  in  kilobars  and  T  is  in  Kelvins.  The  equations  are  valid  for  the  intervals  plot¬ 
ted. 
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this  room-temperature  SWL  curve  is  graphed  as  curve  B  in  figure  1 1. 

3.1.2  HEATED  CELL 

In  order  to  reach  higher  temperatures  behind  the  shock,  the  initial  temperature  of 
the  CS2  sample  was  raised  from  room  temperature  to  i25°C.  A  heated  spectroscopy  cell 
was  devised  for  this  experiment  and  is  shown  in  figure  1 2.  A  heating-coil  groove  was  cut 
around  the  circumference  of  the  brass  body  of  a  standard  cell,  and  a  0.1  inch  o.d.  coil  of 
28  gauge  nichrome  wire  was  cemented  in  the  groove  with  Sauereisen  P-1  refractory 
cement.  This  heater  had  a  resistance  of  2 1  ohms  and  was  powered  with  15-20  volts  AC 
from  a  variable  transformer.  Pre-shot  temperature  was  measured  via  a  thermocouple 
pressed  against  the  cell  body  with  a  screw.  A  washer-shaped  brass  retainer  was  fastened 
over  the  back  window  to  keep  the  press-fit  sapphire  from  being  pushed  out  by  the 
lOOpsia  vapor  pressure  of  CS2  at  125°C.18  In  an  alignment  test  at  125°C,  the  filled  cell 
remained  parallel  to  the  target  ring  face  to  within  0.7  milliradian.  This  value  did  not 
change  after  two  cycles  of  heating  and  cooling.  (Fresh  CS2  was  used  in  the  actual  shock 
experiment) 

A  copy  of  the  original  record  of  the  heated  cell  experiment  is  shown  in  plate  6.  The 
impact  pressure  was  65  kbar  and  the  cell  thickness  was  164  microns.  The  digitized  band 
edge  from  the  original  record  is  graphed  in  figure  1 3  along  with  the  pressure-time  calcula¬ 
tions  from  SHOCK.UP.  AA  -vs-  P  taken  from  this  figure  is  plotted  as  curve  h  in  figure 
14.  The  T-P  profile  for  this  experiment  is  plotted  as  curve  C  in  figure  1 1. 

3.1.3  COOLED  CELL 

A  standard  spectroscopy  cell  was  modified  for  cooling  by  soldering  two  turns  of 
3/16  inch  copper  tubing  around  the  cell  body  prior  to  assembling  the  cell  (figure  15). 
Liquid  nitrogen  was  pumped  through  the  tubing  via  a  long  rubber  hose,  allowing  the  N2 
to  vaporize  before  reaching  the  cell  and  thus  giving  some  measure  of  temperature  control. 
Cell  alignment  was  to  0.4  milliradian  at  room  temperature.  At  -72°C  alignment  was  diffi¬ 
cult  to  measure  due  to  frosting  even  inside  the  glove  bag  built  around  the  alignment  tele- 
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Figure  12:  Heated  Spectroscopy  Ceil 

a)  heater  coil  (insulating  cement  not  shown);  b)  thermocouple  beneath  screw  & 
^Aer’  c  .fhermocouple^leads;  d)  retaining  ring;  e)  fill  tubes;  0  sample  cavity  g)  front 
sapphire  window;  h)  back  sapphire  window;  i)  o-ring  seals;  j)  brass  cell  body. 
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Figure  13:  Heated  Cell  Experiment  (Shot  83-022);  Tq  ■  398 °K 

©  -  AA  -vs-  time;  digitized  band  edge  from  original  shot  record. 

X  -  Pressure  -vs-  time;  calculated  by  SHOCKUP. 

At  time  -  0,  the  shock  just  enters  the  CS5.  At  the  first  X,  the  shock  completes  one  transit 
of  the  CS,  and  is  just  about  to  reflect  from  the  back  sapphire  window.  At  this  instant  the 
CS2  pressure  is  homogeneous  throughout  the  cell.  The  dotted  line  shows  how  pressure, 
calculated  by  SHOCKUP,  is  correlated  with  the  wavelength  shift  of  the  band  edge. 
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Figure  14:  AA  -vs-  P  For  The  3200A  CS-.  Band  Edge 

The  solid  lines  are  least  squares  cubic  fits  to  the  data.  The  equations  of  the  lines  are 
given:  O)  data  from  five  shock  experiments  with  Tn  ■  293°  K,  AA  «  26.75 P  -  0.2554 P2 
+  2.252x  lO^P3.  C)  cooled  cell  experiment,  Tn  -  204°a:,  AA  -  17.404P  -  8.3074x  10‘2P* 
+  1.1 53 1  x  10”3/*3.  h)  heated  cell  experiment/Tn  -  398°*:,  AA  -  45.01 9P  -  0.92222 P3 
+  9.46 15x1  O^P3.  u 

AA  is  in  angstroms,  P  is  in  kilobars. 


Figure  15:  Cooled  Spectroscopy  Cell 

a)  cooling  coil  (3/16th-inch  copper  tubing);  b)  thermocouple  beneath  screw  & 
washer;  c)  thermocouple  leads;  d)  solder;  e)  fill  tubes;  0  sample  cavity;  g)  front  sapphire 
window;  h)  back  sapphire  window;  i)  o-ring  seals;  j)  brass  cell  body. 
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scope,  but  maximum  tilt  was  no  greater  than  ~  I  .0  milliradian. 

Plate  7  shows  a  copy  of  the  original  shot  record  for  the  cooled  cell  experiment  The 
initial  temperature  was  -72°C  (201°K ),  the  impact  pressure  was  61  kbar,  and  the  cell 
thickness  was  158  microns.  The  digitized  band  edge  from  the  original  record  is  graphed 
in  figure  16  along  with  the  pressure-time  calculations  from  SHOCK  UP.  AA  -vs-  P  taken 
from  this  figure  is  plotted  as  curve  C  in  figure  14.  The  T-P  profile  for  this  experiment  is 
plotted  as  curve  A  in  figure  1 1. 

3.2  STATIC  MEASUREMENTS 

In  the  shock  experiments,  spectra  are  recorded  during  compression,  but  also  during 
two  or  three  microseconds  prior  to  impact  Thus  CS2  spectra  were  obtained  at  ambient 
pressure  and  at  20l°K,  293 °K,  and  398°K  (at  398°K  CS2  vapor  pressure  is  6.8  atm., 
which  is  still  ~zero  on  the  scale  of  these  experiments).  See  figure  17.  These  measurements 
show  a  band  edge  shift  of  ~  29 A  as  T  changes  from  201  to  293 °K,  and  ~48Aas  T  changes 
from  293  to  398°AT.  The  spectra  at  ambient  pressure  of  a  100  micron  layer  of  CS2  at 
298°K  and  398 °K  were  also  recorded,  with  a  little  more  accuracy,  on  a  Cary- 14  spectro¬ 
photometer.  In  these  spectra  the  3200 A  band  had  broadened  to  the  red  by  42 A  (at 
Absorbance-2),  as  T  changed  from  29%° K  to  398°K,  but  the  blue  edge  moved  only  6  A  to 
the  blue  (figure  18). 
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Figure  16:  Cooled  Cell  Experiment  (Shot  83-031);  Tq  *  201° AT 

O  -  AA  -vs-  time;  digitized  band  edge  from  original  shot  record. 

X  -  Pressure  -vs-  time;  calculated  by  SHOCKUP. 

At  time  ■  0,  the  shock  just  enters  the  CS*.  At  the  first  X,  the  shock  completes  one  transit 
of  the  CS*  and  is  just  about  to  reflect  from  the  back  sapphire  window.  At  this  instant  the 
CS^  pressure  is  homogeneous  throughout  the  cell.  The  dotted  line  shows  how  pressure, 
calculated  by  SHOCKUP,  is  correlated  with  the  wavelength  shift  of  the  band  edge. 
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Figure  17:  CS2  Spectrum  Showing  Red  Edge  Of  Absorption  Band 

Pre-shock  CS-,  spectra  for  shots  83-031,  82-031,  and  83-022,  at  -72°C  (20 IT), 
20°C  (293T),*and  125°C  (398 °K)  respectively.  These  spectra  were  obtained  by  densi¬ 
tometer  scans  of  the  original  35mm  film  shot  records.  The  figure  shows  the  band-edge 
shift  with  temperature  at  constant  (ambient)  pressure.  The  dotted  line  shows  the  approxi¬ 
mate  level  of  film  density  at  which  all  the  records  were  read.  This  is  the  level  at  which 
the  band  edge  just  becomes  visible  by  eye. 
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Figure  18:  Spectrum  Of  Liquid  CS2 

The  spectra  shown  are  of  a  100  micron  thick  sample  of  CS,  at  298 °K  and  at  398 °K 
recorded  at  ambient  pressure  on  a  Cary  14  spectrophotometer/ 


I 


27 


SECTION  4 

CS2  EQUATION  OF  STATE 

Finding  equations  of  state  for  high  pressures  and  temperatures  has  been  a  long¬ 
standing  problem.  They  are  constructed  from  meager  data,  largely  from  shock  measure¬ 
ments,  which  provide  no  temperature  information.  These  equations  are  commonly 
modelled  after  standard  forms,  for  example,  Woolfolk,  Cowperthwaite,  and  Shaw  have 
proposed  an  empirical  "universal"  liquid  Hugoniotf  based  on  shock  data  for  several 
liquids.  The  CS2  equation  of  state  used  in  this  report  was  written  by  Sheffield  and 
Duvall.20  It  is  based  on  the  Woolfolk,  Cowperthwaite,  and  Shaw  model  which  gives  an 
excellent  fit  to  CS2  shock  data  (pressure  and  volume),  but  gives  no  information  on  tem¬ 
perature.  Sheffield’s  temperature  calculation  is  based  on  an  assumed  function  for  Cy.  To 
give  an  idea  of  the  uncertainties  involved,  the  temperature  of  CS2  after  shock- 
compression  to  50  kbar  has  been  calculated  from  three  different  equations  of  state  (for 
step-wise  loading,  with  impact  pressure  of  54  kbar): 

Sheffield- 1  -  732°K 

Sheffield-II - 610°* 

Mie-Gruneisen - 551  °K\ 

Sheflfield-I  was  written  by  Steve  Sheffield  as  part  of  his  PhD  thesis;  *  *  Sheffield-II  is  the 
Sheffield- Duvall  equation  of  state  used  in  this  report.  It  is  impossible  to  assess  the  accu¬ 
racy  of  these  calculations  without  furthur  thermal  measurements.  The  Sheffield-II  pres- 

21 

sure  calculations  have  been  verified  by  Sutherland,  who  used  a  quartz  gauge  to  meas¬ 
ure  CS2  pressures  up  to  40  kbar  in  an  SWL  (stepwise  loading)  experiment.  The  agreement 
was  within  ~2%. 

For  the  "cold"  experiment  (shot  83-031),  measured  values  of  thermodynamic  state 
functions  for  CS2  at  -72°C  were  obtained  from  the  literature  and  were  used  as  initial 
values  in  the  Sheffield  equation  of  state.  These  quantities  were  V  (specific  volume),  alpha 
(thermal  volume  expansivity),  K-j.  (isothermal  bulk  modulus),  and  Cp  (isobaric  specific 
heat).  The  first  three  were  measured  graphically  from  P-V-T  data  compiled  by  P.  W. 
Bridgman.22  These  data  are  from  different  sources  so  I  will  say  more  about  them.  The 

t  A  Hugoniot  La  the  locua  of  thermodynamic  end-pointa  reached  by  ahock-compreeeion  of  a  material.  It  ia  a 
aingle  curve  on  the  3-dimenaional  thermodynamic  eurface  for  that  material. 

t  WithT/V  conetant,  apecific  heat  waa  allowed  to  vary  with  temperature  according  to  the  Debye  theory. 


•  V  a*  w“  V  Via'1  a" 
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most  complete  P-V-T  measurements  on  CS2  above  0 °C  were  made  by  Bridgman  and 
published  in  1913.  In  these  experiments  T  was  varied  from  0 °C  to  80°C  and  P  was 
varied  from  500  atm  to  12000  atm  (0-500  atm  data  were  taken  from  Amagat,  as  described 
in  Bridgman’s  paper).  These  data  were  extended  below  0 °C  to  -100°C  at  pressures  from 

^  A 

1  atm  to  1000  atm  by  Seitz  and  Lechner  in  1916.  The  Seitz  and  Lechner  data  were 
later  compiled  by  Bridgman  with  his  own  data  and  published  in  the  International  Critical 
Tables.  The  value  for  Cp  at  -72°C  was  taken  from  reference  25.  These  values  are 
listed  in  table  1. 

For  the  125°C  experiment  (shot  83-022),  thermodynamic  quantities  at  ambient  pres¬ 
sure  had  to  be  extrapolated  beyond  known  values,  as  no  measurements  could  be  found 
above  80°C.  V  and  alpha  were  calculated  from  a  cubic  equation  fit  to  V-T  data  for  CS2 
from  -34°C  to  60°C.  This  equation  was  checked  against  the  Bridgman  data  and  the 
Seitz  and  Lechner  data  and  made  an  excellent  fit  from  -34°C  to  the  limit  of  Bridgman’s 
data  at  80°C.  KT  was  calculated  by  plotting  Bridgman’s  data  against  T,  and  extrapolating 
to  125°C.  Cp  at  125°C  was  also  obtained  by  extrapolation,^3  but  should  be  reliable 
because  Cp  changes  little  with  temperature.  The  initial  values  used  in  this  report  are 
listed  in  table  1. 


Table  1 

Values  of  Thermodynamic  Parameters  for  CS2t 


State  Function 

-72°C 

125°C 

20°C 

V  (cm3/g) 

0.7163 

0.9293 

0.7918 

alpha  (°K'1) 

1.0317xl0‘3 

1.9803x1 0'3 

I.2I8xl0"3 

Ky  (Mbar) 

.0211 

.0025 

.0108 

Cp  (Mbar  cmV10*'1) 

0.989x1 0'5 

l.OIxlO*5 

I.OOxIO'5 

t  Thu*  value*  ar*  from  reference*  22  -  26.  Note:  the  125°C  value*  were  extrapolated  beyond  known 
data. 


SECTION  5 

ANALYSIS  OF  RESULTS 


5.1  ANALYSIS 

In  tables  2,  3,  and  4,  I  have  compiled  the  P,  T,  and  AA  data  from  the  heated, 
cooled,  and  room  temperature  experiments  described  in  section  3.  The  values  listed  are 
interpolations  to  convenient  values  of  pressure  and  temperature,  taken  from  the  curves 
for  AA  -vs-  P  and  T  -vs-  P.  (These  interpolations  were  obtained  graphically  and  may 
differ  slightly  from  values  calculated  with  the  equations  fit  to  the  AA  -vs-  P  and  T  -vs-  P 
data,  shown  in  figures  1 1  and  1 4.) 


Table  2 

Cold  SWL  (Shot  83-031) 

Table  3 

Room-Temperature  SWL 

P(kbar) 

AA(X) 

T(°K) 

P(kbar) 

AA(A) 

T(°K) 

(A  *3521  A) 

(A  =3550A) 

0 

0  " 

201 

0 

0  u 

293 

10 

162 

302 

10 

240 

441 

20 

324 

358 

20 

438 

512 

30 

480 

398 

30 

627 

561 

40 

644 

432 

40 

807 

597 

50 

804 

459 

50 

976 

628 

91 - 

~T38 

300 

60 

1142 

656 

30.5 

488 

400 

70 

1373 

679 

46.3 

744 

450 

0.3 

— 7 - 

m~ 

7.1 

175 

400 

18.0 

398 

500 

40.8 

821 

600 

57.7 

1095 

650 

70.5 

1384 

680 

From  the  upper  part  of  each  table,  AX/ AT  can  be  calculated  at  convenient  pressures,  and 
from  the  lower  part  AX/A P  is  obtained  at  convenient  temperatures.  For  example, 
AX /AT  at  10  kbar,  between  cold  and  room-temperature  experiments: 

Cold  experiment :  A  -  AA  +  XQ  =  162A  +  3521  A  =-  3683 A 
Room-temp.  exp. :  X  -  AA  +  XQ  -  24oA  +  3550A  -  3790A 


AA  _  3790  -  3683  _  107A 

AF  J10  *  441  -  302  139°A: 


0.77  A/°K 


Thus,  we  can  construct  tables  5  &  6.  The  data  from  these  tables  are  plotted  in  figures  19 
through  22. 


Figure  19:  AA/AT  -vs-  T 

X  -  static  data  (sec.  3.2) 

+  -  between  cold  and  room-temperature  experiments 
X  -  between  cold  and  hot  experiments 
O  -  between  room-temperature  and  hot  experiments 
The  line  is  a  least  squares  fit  of  all  the  data  and  has  the  equation: 

~  -  0.415  +  0.003747(r-300°/i:) 

where  AA/AT  =  k/°K. 


Figure  20:  AA  /  AT  -vs-  P 

X  -  static  data  (sec.  3.2) 

+  -  between  cold  and  room-temperature  experiments 
X  -  between  cold  and  hot  experiments 
O  -  between  room-temperature  and  hot  experiments 
The  number  next  to  each  point  is  the  average  temperature  in  °K  for  that  point  (from  table 
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Figure  22:  AA/AP  -vs-  P 

+  -  between  cold  and  room-temperature  experiments 
X  -  between  cold  and  hot  experiments 
O  -  between  room-temperature  and  hot  experiments 


dX/dl  and  ax/af : 


If  we  assume  that  dX/dT  is  given  by  the  line  drawn  in  figure  19,  then  we  can  calcu¬ 
late  dX/dP  for  each  experiment  as  follows: 


d  X 


ax 

ar 


dT  + 


ax 

ap 


dP  ( assume  X  -  X(T,P)) 


ax  ' 

dX  | 

ax 

ap 

T  "  dP  “  \ 

dT 

dT 

pdP 


For  example,  AA  -vs-  P  for  shot  83-031  (cold  SWL)  is  plotted  as  curve  C  in  figure  14. 
The  tangent  at  any  point  along  the  curve  gives  dX/dP.  From  figure  19,  dX /dT  =*  0.4 15  + 
0.003747(T  -  30(FK ),  and  dT/dP  is  the  slope  of  the  tangent  to  the  T  -vs-  P  curve  in  fig¬ 
ure  11.  Thus  we  can  calculate  dX/dP  for  any  pressure  in  an  experiment  The  values  of 
dX  /dP  calculated  at  10  kbar  intervals  for  the  hot,  cold,  and  room-temperature  experi¬ 
ments  are  plotted  against  temperature  and  pressure  in  figures  23  and  24  respectively. 
dX/dP  displays  no  dependence  on  either  temperature  or  pressure.  The  average  value  for 
dX/dP  is  13.4  A/ kbar. 


Self-Consistency  Test : 

The  above  values  for  dX  /dT  and  dX/dP  are  applied  to  each  experiment  to  try  and 
reproduce  the  AA  -vs-  P  data  by  using  the  following  relationship: 

AA  -  |^A T  +  |^A P 
dT  dP 

The  results  of  this  test  are  listed  in  table  7.  A  slightly  better  overall  fit  is  obtained  with  a 
value  of  12.8  kjkbar  for  dX/dP. 
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Figure  23:  3X/3P  -vs-  T 

+  -  cooled  cell  experiment  (shot  83-031) 

*  -  room  temperature  experiment  (shot  82-037) 

©  -  heated  cell  experiment  (shot  83-022) 

dX/d?  was  calculated  from:  —  »  —  -  Q-dJ. 

dP  dP  dT  dP 

where i  dX/dT  -  0.415  +  ,003747(T  -  200°K)  (from  fig.  19).  and  dT/dP  and  dA/dP  are 
trom  figures  1 1  and  14  respectively.  A  is  in  angstroms,  P  is  in  kilobars,  and  T  is  in  Kel- 
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Figure  24:  3A/3P  -vs-  P 

+  -  cooled  cell  experiment  (shot  83-031) 

di  -  room  temperature  experiment  (shot  82-037) 

O  -  heated  cell  experiment  (shot  83-022) 

_.  /an  -  -  „  .  -  dX  d X  dX  dT 

dX/dP  was  calculated  from:  —  *  ~zzr  “  — — 

dP  dP  dT  dP 

where  dX/dT  -  0.415  +  ,003747(T  -  300oA:)  (from  fig.  19),  and  dT/dP  and  dA/dP  are 
from  figures  11  and  14  respectively.  A  is  in  angstroms,  P  is  in  kilobars,  and  T  is  in  Kel- 


Tabic  7  :  Seif-Consistency  Test:  Comparison  Of  Measured  Values  of  AA  (Tables  2-4)  With 
Calculated  Values  From  The  Equation: 


AA 


3\a  t  +  44a  P 


3T  3P 

with  3X/3T  -  0.415  +  ,003747(T  -  300°AO,  and  dX/dP  -  13.4  (A  is  in  X,  T  is  in  °K,  P  is  in 
kbars). 


Cold  SWL  (SHOT  83-031) 


AA 


A  P(kbar) 

A  T(°K) 

dX/dP 
( Jt/kbar ) 

dX/dT 

(A/°K) 

calculated 

(A) 

actual 

(A) 

error 

% 

0-10 

201-302 

13.4 

0.23 

157.3 

162 

-2.9 

10-20 

302-358 

13.4 

0.53 

163.5 

162 

+0.9 

20-30 

358-398 

13.4 

0.71 

162.3 

156 

+4.0 

30-40 

398-432 

13.4 

0.85 

162.8 

164 

+0.7 

40-50 

432-459 

13.4 

0.96 

160.0 

160 

0.0 

Room-Temperature  SWL 

0-10 

293-441 

13.4 

0.67 

232.6 

240 

-3.1 

10-20 

441-512 

13.4 

1.07 

210.3 

198 

+6.2 

20-30 

512-561 

13.4 

1.30 

197.6 

189 

+4.6 

30-40 

561-597 

13.4 

1.46 

186.6 

180 

+3.6 

40-50 

597-628 

13.4 

1.58 

183.1 

169 

+8.3 

50-60 

628-656 

13.4 

1.70 

181.5 

166 

+9.3 

Hot  SWL  (SHOT  83-022) 

0-10 

398-586 

13.4 

1.13 

347.2 

345 

+0.6 

10-20 

586-644 

13.4 

1.60 

226.5 

258 

-12.2 

20-30 

644-681 

13.4 

1.77 

199.5 

187 

+6.7 

30-40 

681-709 

13.4 

1.90 

187.1 

161 

+16.2 

40-50 

709-731 

13.4 

1.99 

177.8 

170 

+4.6 

5.2  UNCERTAINTIES 


The  pressures  and  temperatures  calculated  in  this  report  are  based  on  the 
ShefTield-II  equation  of  state.  The  accuracy  of  this  equation  to  within  2%  for  shock- 
pressure  calculations  has  been  verified  up  to  40  kbar.21  Temperature  calculations,  how¬ 
ever,  depend  on  an  assumed  function  for  Cy  with  shock  pressure.  Different  equations  of 
state  can  give  widely  differing  results.  For  example,  Lysne  and  Hardesty27  have  shown 
that  calculated  shock  temperatures  for  nitromethane  at  200kbar,  can  differ  by  nearly  a 
factor  of  two  for  different  equations  of  state.  The  two  Sheffield  forms  (section  4)  predict 
temperatures  of  732 °K  and  610°A:  at  50kbar  in  a  SWL  experiment  The  later  Sheffield 
form,  given  in  reference  20,  was  used  in  this  report 

There  may  be  unknown  errors  introduced  into  the  Sheffield  equation  of  state  by 
changing  Tq.  Other  uncertainties  arise  in  producing  and  reading  the  photographic  shot 
records.  The  reproducibility  in  obtaining  the  pre-shock  position  of  the  CS2  band  edge  in 
repeat  experiments  is  ~±15/£.  In  correlating  the  wavelength  shift  with  the  calculated 
shock  pressure,  as  in  figures  9,  13,  and  16,  small  errors  in  ceil  thickness  measurement  can 
give  rise  to  significant  errors  in  the  calculated  time-position  of  the  reverberating  shock 


wave. 


SECTION  6 

OTHER  MEASUREMENTS 
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6.1  DIAMOND  CELL 

The  3200A  absorption  band  in  CS2  has  been  studied  in  a  diamond  anvil  cell  by 
28 

Agnew  &  Swanson.  They  measured  the  position  of  the  red  edge  of  this  band  as  a  func¬ 
tion  of  static  pressure  at  ambient  temperature.  Their  results  are  plotted  in  figure  25  along 
with  the  shock  measurements  from  this  report  For  the  diamond  cell  data,  dX/dP  ■  9.9 
kjkbar ,  which  is  26%  less  than  the  shock  data  average  of  13.4  k/kbar. 

6.2  NOTE  ON  REACTION  THRESHOLD 

Dick  ’  u  found  a  transition  to  a  higher  density  state  in  CS2  occurring  in  single¬ 
shock  experiments  at  about  60  kbar.  The  measurements  made  in  the  present  report  are  on 
"unreacted"  CS2,  that  is,  below  the  temperature  and  pressure  threshold  of  the  transition 

observed  by  Dick.  In  stepwise  loading  experiments,  Duvall,  Bellamy,  Hegland,  Eschbach, 

29 

and  this  author  showed  that  an  irreversible  spectral  change  occurs  in  CS2  at  89  kbar. 
This  change  appears  to  be  the  same  transition  observed  by  Dick,  but  occurs  at  a  higher 
pressure  due  to  the  lower  temperature  in  a  multiple-shock  compression  as  opposed  to  a 
single-shock  compression. 


Figure  25:  P-T-A  Surface  For  Liquid  CS2 

This  surface  shows  the  position  of  the  red  edge  of  the  3200A  absorption  band  in  CS?,  as 
a  function  of  pressure  and  temperature.  Three  of  the  curves  depicting  the  surface  are-  the 
results  of  shock  experiments  on  CS^at  three  different  initial  temperatures:  201  °K,  293 °K, 
and  398 °K.  The  isothermal  curve  at  293 °K  is  the  result  of  diamond  anvil  cell  measure¬ 
ments  made  by  Agnew  and  Swanson.38  The  dotted  lines  connect  points  of  constant  pres¬ 
sure. 


SECTION  7 

CONCLUSIONS  /  RECOMMENDATIONS 

Heated  and  cooled  spectroscopy  cells  were  successfully  employed  in  shock  experi¬ 
ments  to  recoru  .  ^e-resolved  spectra  of  liquid  CS2.  The  red  edge  of  the  lowest  energy 
electronic  absorption  band  in  liquid  CS2  exhibits  a  highly  regular  shift  that  appears  to  be 
a  function  only  of  pressure  and  temperature.  The  temperature  component  seems  indepen¬ 
dent  of  pressure,  but  appears  to  increase  linearly  with  temperature  by  dX/dT  *  0.415  + 
0.003747(T  -  300°K).  The  pressure  component  of  the  shift  appears  roughly  constant  at 
~13  k/kbar  ±~20%.  Static  data  (diamond  cell)^  indicate  a  value  of  9.9  A/kbar  for 
isothermal  compression  at  room  temperature.  Further  thermal  measurements  need  to  be 
made  on  CS2  at  high  pressures  and  temperatures  in  order  to  remove  the  uncertainty  in 
the  equation  of  state  temperature  calculations,  and  to  provide  coirect  thermodynamic 
values  for  initial  conditions  in  the  experiment  at  125°C. 
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APPENDIX 


SUMMARY  OF  EXPERIMENTS 


Shot  Number 

T0(°K) 

Impact  Pressure 
(kbar) 

Final  T(°K) 

Cell  Thickness 
(microns) 

82-031 

293 

54 

617 

165 

82-028 

293 

75 

683 

116 

82-037 

293 

81 

701 

160 

83-014 

293 

130 

844 

158 

82-015 

293 

71 

670 

180 

83-022 

398 

65 

757 

164 

83-031 

201 

61 

486 

158 

A 


y>; 


,  -*. 

■  ^  — 


PLATE  1 


|  1  MS 

impact 


SHOT  82-031 

Impact  pressure  a  54  kbar 
Cell  thickness  a  165  microns 


PLATE  2 


impact 


SHOT  82-028 

Impact  pressure  »  75  kbar 
Ceil  thickness  s  116  microns 


impact 


SHOT  82-037 

Impact  pressure  a  31  kbar 
Cell  thickness  ■  160  microns 
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PLATE  4 


Time 


f  1  MS 


impact 


SHOT  83-014 

Impact  pressure  a  130  kbar 
Cell  thickness  a  158  microns 


(Apparent  initial  absorption  band  edge  is  at  3710A  This  is  higher  than  the  true  band 
edge  of  3550a  becq  o  in  this  experiment  the  light  passed  through  a  lucite  filter  which 
absorbs  up  to  3710A 


PLATE  5 


SHOT  82-015 


Impact  pressure  *  71  kbar 
Cell  thickness  ■  180  microns 

Note:  the  wavelength  scale  of  this  record  is  expanded  4X. 


impact 


SHOT  83-031 


Cooled  cell  experiment 
Impact  pressure  ■  61  kbar 
Cell  thickness  *  158  microns 
>  201°K 
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